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Abstract 

For a simple and connected graph, a new graph invariant s* (G), defined as 
the sum of a powers of the eigenvalues of the normalized Laplacian matrix, has 
been introduced by Bozkurt and Bozkurt in [7j. Lower and upper bounds have 
been proposed by the authors. In this paper, we localize the eigenvalues of the 
normalized Laplacian matrix by adapting a theoretical method, proposed in Bianchi 
and Torriero ( 0 ), based on majorization techniques. Through this approach we 
derive upper and lower bounds of s*(G). Some numerical examples show how 
sharper results can be obtained with respect to those existing in literature. 


1 Introduction 

Among the various indices in Mathematical Chemistry, a whole new family of descriptors 
s*(G), defined as the sum of a powers of the eigenvalues of the normalized Laplacian 
matrix, has been proposed by Bozkurt and Bozkurt in [£j. These authors found a number 
of bounds for arbitrary a and particularly for a = — 1, which is the case of the degree 
Kirchhoff Index. Recently, Bianchi et al. proposed a variety of lower and upper bounds 
for s* a (G) in |lj and for the Kirchhoff Index in [2j derived via majorization techniques. 
In particular, the authors showed that it is possible to obtain tighter results taking into 
account additional information on the localization of the eigenvalues of proper matrices 
associated to the graph. From a theoretical point of view, some well-known inequalities 


on the localization of real eigenvalues have been provided in literature and they can be 
used to compute the above mentioned bounds. 

Alternative inequalities involving the localization of some eigenvalues of the transition ma¬ 
trix of the graph have been numerically computed in [9j and [TO} by adapting a theoretical 
methodology proposed in Bianchi and Torriero [5] based on nonlinear global optimization 
problems solved through majorization techniques. By means of these results, tighter lower 
bounds for the Kirchhoff Index for some classes of graphs have been derived in p]. 

The purpose of this paper is to use this fruitful theoretical method providing indeed some 
formulae that allow us to compute lower bounds for the first and the second eigenvalues 
of the normalized Laplacian matrix in a fairly straightforward way. Thus, we obtain new 
bounds for s* (G) considering both non-bipartite and bipartite graphs. 

In Section 2 some preliminaries are given. In Section 3 we describe the nonlinear 
optimization problem useful for our analysis. Solving this optimization problem, lower 
bounds of the first and second eigenvalue of the normalized Laplacian matrix have been 
provided in Section 4. Finally, in Section 5 several numerical results are reported showing 
how the proposed bounds are tighter than those given in [7]. 

2 Notations and preliminaries 

In this section we first recall some basic notions on graph theory. For more details refer 
to [17] . 

Let G = (V, E) be a simple, connected, undirected graph where V = {1, 2,..., n} is the 
set of vertices and E C V x V the set of edges, \E\ = m. 

The degree sequence of G is denoted by tt = (d 1; d 2 , d n ) and it is arranged in non¬ 
increasing order d\ > d 2 > • • • > d n , where dj is the degree of vertex i. It is well 

n 

known that = 2m and that if G is a tree, i.e. a connected graph without cycles, 

i=l 

m — n — 1. Let A(G) be the adjacency matrix of G and D{G) be the diagonal matrix 
of vertex degrees. The matrix L(G) = D(G ) — A(G) is called Laplacian matrix of G, 
while E(G) = D(G)- 1 / 2 L(G)D(G)~ 1 / 2 is known as normalized Laplacian matrix. Let 
hi ^ ^2 > > hn be the (real) eigenvalues of L{G) and Ai > A 2 > ... > A n be the (real) 

eigenvalues of C(G). The following properties of spectra of L{G) and C(G) hold: 

n ' 2m 

/Pi = tr (L(G)) = 2m; hi > 1 + d i > —; hn = 0, Hn-i > 0. 
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It Ifc 1 

£> = tr(£(G)) = n; £\\f = tr(£ 2 (G)) = n + 2 £ —; A„ = 0;A,<2. 

»=1 *=1 ( ij ) e £ * 3 

Our aim is the analysis of a particular topological index, s* (G). Topological indices have 
been widely explored in different fields, i.e. mathematical chemistry and more recently in 
complex network analysis. In particular, Zhou (see [18]) proposed the index: 

71—1 

s a {G) = ± 0,1, 

2—1 

defined as the sum of the a-th power of the non-zero Laplacian eigenvalues of a graph G. 

Over the last years this index and its bounds have been intensely studied: Zhou (see 
[T8] ) established some properties of s a (G) and some improvements have been provided in 
[H], [16], [19] and [20]. In [3], taking into account the Schur-convexity or Schur-concavity 
of the functions s a (G ) for a > 1 and a<0or0<a<l respectively, the same bounds 
as in [18] have been derived. Furthermore, considering additional information on the 
localization of the eigenvalues, the authors provide also sharper bounds. 

Bozkurt and Bozkurt in [7] introduced parallely to [18] the following new graph in¬ 
variant: 

n— 1 

4(G) =£v, a/0,1, 

2=1 

characterized as the sum of the a-th power of the non-zero normalized Laplacian eigen¬ 
values of a graph. Several properties of this index have been proposed in [7] and some 
lower and upper bounds for a connected graph have been derived. 

In [1], considering the Schur-convexity or Schur-concavity of the functions s*(G) and 
using additional information on the localization of the eigenvalues, the following Theo¬ 
rems, which generalize Theorem 3.3 in [7], have been proved. 

Theorem 1. Let G be a simple connected graph with n > 3 vertices and X\ >9: 

1. if a < 0 or a > 1 then 

S “ (G)S9 " + (fe>FT « 

2. if 0 < a < 1 then 

+ < 2 > 

Theorem 2. Let G be a simple connected graph with n > 4 vertices which is not complete 
and Ai > 6, \ 2 > /3 with 6 > f3 and 9 + j3(n — 2) > n. 
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1. if a < 0 or a > 1 then 


s* a (G) >O a + p a + 


(n — 9 — /3) a 
(■n - 3 ) Q_1 


(3) 


2. if 0 < a < 1 then 

4(G) <r + fS‘+ (4) 

It is noteworthy to state that the results in Theorem [2] are tighter than those in 
Theorem [I] (for more details see [3] and W)- 

In [7], the bounds in Theorem [I] has been previously proved identifying 9 as 


P = 


1 + 


n\n 


1) 


E 

hhe-B) 


1 

didj 


In Section [ij by adapting some results proved in |5], we provide lower bounds of Ai 
and A 2 that enable us to obtain tighter bounds than in [7]. In what follows we refer to the 
values of lower bounds of Ai and A 2 as Q and R , respectively, i.e. Ai > Q and A 2 > R. 


3 A Nonlinear optimization problem to bound eigen¬ 
values 


We now recall a methodology based on majorization techniques (see [5] and ra) that 
allows us to find a suitable localization of Ai and A 2 in order to provide some new results 
for bounds of s* (G). 

At this regard, we define the set 

n—1 n—1 

S* = {A e M^ 1 : Ar > A 2 > ... > X n -i > 0, ^ A, = n,g( A) = ^Af = b}, 

2—1 2=1 

where p is an integer greater than 1 . 

The following fundamental lemma holds (see Lemma 2.1 in [IJ]): 


Lemma 1. Fix b e M and consider the set Sf. Then either b 
unique integer 1 < h* < (n — 1 ) such that: 


— 7 —or there exists a 

(n— 1)P 1 


n F 


(h* + 1 ) 


^I< b < 


n p 


( h*) 


p— 1 5 


(5) 


where h* = p ~f X- 
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The following Theorem is derived by Theorem 3.2 in [5j. Moreover, a lower bound for 
A h by solving the following optimization problem: 


min ( X h ) subject to A 6 5^ 


P*(h) 


Theorem 3. The solution of the optimization problem P*(h) is (^-) if b — ^ n Jiy P -i 


If b ^ (yzyywr. the solution of the optimization problem P*(h ) is 5* where 

1. for h — 1, S* is the unique root of the equation 

f(5, p) = h*5 p + (n- h*5) p -6 = 0 

2. for 1 < h < (h* + 1), 5* is the unique root of the equation 


( 6 ) 


f{6,p) = {n - h)5 p + (h- 1) ^ — 6 = 0 


(h-iy 


ml = ( 0 ,^]; 

3. for h> (h* + 1), 5* is zero. 


(7) 


4 New bounds for normalized Laplacian eigenvalues 


In order to find new bounds for s* (G), we make use of the methodology introduced in 
Section [3] that allows us to provide lower bounds for Ai and A 2 . At this regard, we now 
consider Theorem [3] limiting the analysis when p — 2: in this case we know indeed that 

6 = n + 2 Yl(i,j)eE d~dj- 
2 

When 6 = then the solution of optimization problem is ^zrp- This is the case 

2 2 

of the complete graph K n . Instead, when 6 ^ h* — \ ■ 

Since we want to get a lower bound for Ai, we solve equation <§- The acceptable 

solution in the interval I is equal to 5* = —- n+h*) - - an h we re f er to this value as Q. 

The value of Q can be compared to P in order to show how bounds ([Tj) and ([2]), 
computed by assuming 6 = Q, perform better than those for 9 = P. 

It is well known that, for every connected graph of order n (see 0), we have: 


1 

n — 1 



n 


E 


1 

di dj 


< 1, 


( 8 ) 
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and the left inequality is attained for the complete graph G = K n . 

Figure jlj reports patterns of P and Q , varying the quantity t — 2 • ge dir ^ ie 
proper interval (^-, 71 ) (see Equation (|8 )) for alternative values of number of vertices n. 
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Figure 1: Q and P according to different values of t G (^-p,n) and several number of 
vertices. 

Being P = 1 + y, it is easy to see that P has a monotonic behaviour with 

respect to t in the interval 1 + \J yzy) • Figure jlj shows that Q assumes values equal 

to P when t = (indeed P = Q = in this case) and then it increases faster than 
P in the considered interval. 

Going deeply into the analysis, we graphically show that the difference Q — P is always 
strictly greater than zero according to several values of n and t (see Figure |2|. 

Since Ai > Q > P, we provide bounds (JTj) and ([2]) better than in [7j (see [3) and |!j 
for more theoretical details). 

With the aim to improve previous results, we can now derive additional information 
on A 2 . To evaluate the lower bound R of A 2 , from Theorem |3j considering the case h = 2, 
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since 


h < (h* + 1), we solve the equation ([7]) finding the acceptable solution: 


n — 


R = 5* = 


b(n—l)—n 2 


n —2 


n — l 


In this case, the rightmost inequality in ([8]) implies t < n and then b < 2 n. By 
plugging this information in the value of R , we easily obtain R < that fulfills the 
condition R < Q of Theorem [2j The other condition of Theorem [2] , Q + R(n — 2 ) > n, 
will be numerically checked in the next Section. In what follows, fixed values of Q and 
R, we can numerically compute bounds ([3]) and Q and comparing with those in [7J. 

In case of bipartite graphs it is well known that Ai = 2. If we set 6 = 2 in Theorem 
[lj we derive the same results found in [Ij. Furthermore, by placing 6 = 2 and (3 = R in 
bounds ([3]) and Q, we provide bounds also for bipartite graphs. 


5 Some Numerical Results 


The proposed bounds have been evaluated on different graphs. We now focus only on 
non-bipartite graphs and we provide a comparison with literature (see 0). 

In order to assure a robust analysis, graphs have been randomly generated following the 
Erdos-Renyi (ER) model GEn(’n,q ) (see [6j, [8], [llj and p2j). Graphs have been obtained 
by using a MatLab code that gives back only connected graph based on the ER model 
(see |9j and [TO]). In this fashion, the graph is constructed by connecting nodes randomly 
such that edges are included with probability independent from every other edge. The 
results are based on a classic assumption of a probability of existence of edges q equal to 
0.5. We obtain indeed that the generated graphs have a number of edges not far from the 
half of its maximum value as proved in the literature (see for example m)- 

At this regard, in Table [lj s* (G) has been computed for several graphs by fixing 


7 









or 


a equal to 0.5. We report values of upper bound Q evaluated by using 9 = Q 
6 = P (as proposed in ®)- We refer to these bounds as (|2p) and (|2]P). Likewise bound 
(|4}5-R) identifies bound Q evaluated when 9 = Q and (3 = R, where the results has been 
provided assuring that assumptions of Theorem [2] are satisfied. Relative errors r measures 
the absolute value of the difference between the upper bounds and s* a (G) divided by the 
value of s* (G). We observe an improvement with respect to existing bounds according 
to all the analyzed graphs and the improvement appears reduced for very large graphs. 
However, for large graphs the formula provided in [7] already gives a very low relative 
error. 


n 

dl 

m 

s*(G) 

bound |2^>) 

bound |4[3-R) 

bound |2jP) 

r|2]2) 

r|4]3H) 

r|2|p) 

4 

2 

3 

3.35 

3.44 

3.43 

3.46 

2.86% 

2.55% 

3.47% 

5 

4 

9 

4.46 

4.47 

4.47 

4.47 

0.23% 

0.21% 

0.25% 

6 

3 

6 

5.30 

5.47 

5.46 

5.48 

3.13% 

3.00% 

3.27% 

7 

5 

14 

6.43 

6.48 

6.48 

6.48 

0.83% 

0.81% 

0.86% 

8 

5 

13 

7.33 

7.48 

7.48 

7.48 

2.02% 

1.98% 

2.06% 

9 

6 

16 

8.31 

8.48 

8.48 

8.48 

2.04% 

2.01% 

2.07% 

10 

8 

25 

9.39 

9.51 

9.48 

9.52 

1.36% 

1.04% 

1.37% 

20 

15 

95 

19.37 

19.51 

19.49 

19.51 

0.71% 

0.62% 

0.72% 

30 

19 

209 

29.36 

29.50 

29.50 

29.50 

0.49% 

0.46% 

0.49% 

50 

33 

604 

49.37 

49.50 

49.50 

49.50 

0.27% 

0.26% 

0.27% 

100 

60 

2459 

99.37 

99.50 

99.50 

99.50 

0.13% 

0.12% 

0.13% 

200 

116 

10001 

199.38 

199.50 

199.50 

199.50 

0.06% 

0.05% 

0.06% 

300 

179 

22437 

299.37 

299.50 

299.50 

299.50 

0.04% 

0.04% 

0.04% 

500 

279 

62456 

499.38 

499.50 

499.50 

499.50 

0.03% 

0.02% 

0.03% 


Table 1: Upper bounds for s* (G) for a = 0.5 and relative errors. 


The comparison has been extended in order to test the behaviour of the upper bounds 
on alternative graphs. First of all, in the ER model used to generate graphs, the parameter 
q can be thought of as a weighting function. As q increases from 0 to 1, the model becomes 
more and more likely to include graphs with more edges and less and less likely to include 
graphs with fewer edges. In this regard, we assign several values of q moving from the 
default value of 0.5. For sake of simplicity we report only the relative errors derived for 
graphs generated by using respectively q = 0.1 and q = 0.9 (see Figure [3]). In all cases 
bound (|4}5R) assures the best approximation to s* a (G ) for a = 0.5. We observe a best 
behaviour of all bounds when q = 0.9 because we are moving towards the complete graph. 
We have indeed that the density of the graphs increases as long as greater probabilities 
are considered. 

Finally, for the same index So. 5 (G), upper bounds have been evaluated for tree^j] Table 
[2] depicts slighter differences for larger graphs in this case too. However it could be noticed 

1 Tree has been generated by using Priifer code. The Prater sequence of a labeled tree is a unique 
sequence associated to the tree. The sequence for a tree on n vertices has length n — 2 and it can be 
generated by a simple iterative algorithm. It is a way to map bijectively trees on n vertices into n — 2 
long sequences of integers drawn from n. 
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Figure 3: Relative errors of upper bounds of Sq 5 (G) for graphs ER(n, 0.1) and ER{n , 0.9) 
respectively. 


how the relative improvement of bounds respect to other bounds is greater than in case 
of non-bipartite graphs. Despite greater relative errors are observed, bound ([dJ^R) is 
confirmed as the tighter bound also in this case. 



Trees 

n 

<C(G) 

rfSp) 

r|4pfl) 

r|2f) 

4 

3.35 

2.04% 

1.95% 

3.47% 

5 

4.32 

2.17% 

2.14% 

3.48% 

6 

5.23 

3.56% 

3.51% 

4.73% 

7 

6.19 

3.62% 

3.59% 

4.67% 

8 

7.15 

3.67% 

3.64% 

4.61% 

9 

8.22 

2.35% 

2.34% 

3.20% 

10 

8.85 

6.35% 

6.32% 

7.15% 

20 

18.07 

7.45% 

7.44% 

7.88% 

30 

27.63 

6.47% 

6.47% 

6.77% 

50 

45.73 

8.07% 

8.06% 

8.25% 

100 

91.23 

8.97% 

8.97% 

9.06% 

200 

182.72 

9.14% 

9.14% 

9.18% 

300 

274.71 

8.99% 

8.99% 

9.02% 

500 

457.71 

9.11% 

9.11% 

9.13% 


Table 2: So 5 (G) and relative errors for Trees T. 


The analysis has been further developed considering a value of a equal to 1.5. Gener¬ 
ating a similar sample of graphs, both 5 (G) and the relative bounds have been derived. 
For sake of simplicity we report only the results for ER(n, 0.5) observing that the addi¬ 
tional information on the localization of Ai and A 2 lead to the tighter lower bound (|3]5-R). 
Analogous results have been obtained by considering both ER graphs with alternative 
values of q and bipartite graphs. 
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n 

di 

m 

<(G) 

bound |^3) 

bound |3|3-R) 

bound jl^ 1 ) 

r 0n 


r{T]p) 

4 

3 

4 

4.79 

4.66 

4.67 

4.62 

2.69% 

2.39% 

3.49% 

5 

2 

4 

6.22 

5.65 

5.69 

5.60 

9.07% 

8.51% 

9.95% 

6 

4 

9 

6.85 

6.59 

6.60 

6.57 

3.78% 

3.63% 

4.00% 

7 

6 

13 

7.77 

7.57 

7.57 

7.56 

2.56% 

2.49% 

2.65% 

8 

7 

18 

8.75 

8.56 

8.56 

8.55 

2.15% 

2.10% 

2.20% 

9 

4 

12 

10.28 

9.58 

9.59 

9.55 

6.88% 

6.79% 

7.15% 

10 

7 

26 

10.83 

10.52 

10.55 

10.51 

2.90% 

2.60% 

2.94% 

20 

13 

98 

20.88 

20.50 

20.52 

20.50 

1.81% 

1.71% 

1.81% 

30 

19 

222 

30.88 

30.50 

30.51 

30.50 

1.21% 

1.18% 

1.21% 

50 

31 

644 

50.85 

50.50 

50.51 

50.50 

0.68% 

0.67% 

0.68% 

100 

62 

2512 

100.87 

100.50 

100.50 

100.50 

0.37% 

0.36% 

0.37% 

200 

117 

9918 

200.88 

200.50 

200.50 

200.50 

0.19% 

0.19% 

0.19% 

300 

179 

22540 

300.87 

300.50 

300.50 

300.50 

0.12% 

0.12% 

0.12% 

500 

279 

62063 

500.88 

500.50 

500.50 

500.50 

0.08% 

0.08% 

0.08% 


Table 3: Lower bounds for sj 5 (G) and absolute value of relative errors. 


6 Conclusions 


In this paper we provide tighter bounds for the sum of the a-power of the non-zero normal¬ 
ized Laplacian eigenvalues taking into account additional information on the localization 
of the eigenvalues of the normalized Laplacian matrix of the graph, C(G). To this aim 
lower bounds of the eigenvalues are derived by means of the solution of a class of suitable 
nonlinear optimization problems based on majorization techniques. We provide indeed 
closed formulae that allow to compute upper and lower bounds of s* (G) by using the 
additional information on the first and the second eigenvalue of C{G). Numerical com¬ 
parisons confirm how bounds are former than those existing in literature. In particular, 
the analysis has been developed randomly generating both bipartite and non-bipartite 
graphs with a different number of vertices. 
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